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Abstract 

We present a tight-binding calculation of the influence of distorsion on the bulk electronic struc- 
ture of the chalcopyrite CuInSe2- We calculate the ideal case and then the effect of the inclusion 
of the distortions. We analyze our results in detail and conclude from a comparison with other 
work that the distortions must be included in the Hamiltonian to get a proper account of the elec- 
tronic band structure. We use our new Hamiltonian to study the effect that both the tetragonal 
and the anionic distortion have on the (112) surface electronic band structure. We find this effect 
non- negligible.. 
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I. INTRODUCTION 



The quest for room temperature ferromagnetic semiconductors that can be matched to 
conventional semiconductors resulted in an increasing interest in A 7i B /y C^ as well as in 
A J B m C^ chalcopyrites [J. These materials are also interesting as non-lineal optic devices, 
infraredphotodetectors and in solar cell applications with a high efficiency-to-cost ratio [2, 

DDQDE 

Chalcopyrites are tetragonal centered crystallographic structures with eight atoms in the 
unit cell basis. Their spatial group is the D^j. The location and identification of the eight 
atoms in the CuInSe2 unit cell is shown in Table IV below. 

Chalcopyrites can deviate from their ideal symmetry in two ways. The c/a ratio can be 
different from its ideal value 2 and also the anion which lies in the middle of a tetrahedron 
can slide along the central axis. In this paper, we present a detailed study of a tight-binding 
calculation that includes the influence that both distortions have on the bulk electronic band 
structure of CuInSe2 . It is straightforward to generalize this work to calculate their influence 
in the whole series of Cu-based chalcopyrites. The final result is a set of Hamiltonians precise 
enough for the whole series that allow to study different surfaces, monolayers and interfaces 
of these materials in more complicated systems and situations of current interest with a very 
simple method. 



II. SOME DETAILS ON THE METHOD 



re distorted 



10, 
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We use the tight-binding method [8] to calculate the non-distorted as well as t 
case. This method has been used before to describe chalcopyrites successfully j^, 
as well as the related zincblende semiconductor compounds 14 1. A / B /// C 2 // 

In the tight-binding method one constructs Bloch functions 0^(k, r) to describe an elec- 
tronic orbital v centered at the position r+d^ of the ion labelled by fi, as a linear combination 
of atomic-like orbitals ipv( r ) 



J£(k,r) 



/ a z - 



k.(r+d M )./,^ 



<(r-(r + d,)) 



(1) 



where k is a Bloch vector in the First Brillouin Zone and N the number of unit cells in the 
crystal volume considered. We describe the group III metal In and the group IV anion with 
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atomic-like orbitals of s- and p 3 — symmetry. For Cu we consider a full s, p 3 , d 5 basis. The 
spin-orbit interaction is responsible for a crystal field splitting of the heavy and light hole 
bands on the top of the valence band. It is not taken into account in our calculation since 
it is expected to be small. The matrix elements of the Hamiltonian have the form: 



J ^(k, v)H^,(k', r')dr = 6 K « £ e ik ' d -' < v \ v > ><w (2) 

d,, ,,/ 



where 



< v\v >d 



(3) 



dun' is the position vector of the // atom from the a atom. To calculate the non-diagonal 



matrix elements in El we use the Harrison's rule 15[. Therefore the interaction between 
an atomic-like orbital of symmetry x located at the site \x = 1 (In) with another atom of 
symmetry y at // = 2 (Se) is given by = lm[V (ppa) — V (ppn)] . To actually calculate the 
tight-binding parameters, we use further V(ija) = r](ija)H 2 /md 2 l /(d^/ is the interatomic 
distance, m the electron bare mass) for s and p atomic-like orbitals. For the interaction 



between s, p 
III and ref. 



with d orbitals, we use instead V(ida) = rj(ida) h 2 r^ 2 / md'Jf, (for see Table 



j7/ 2 



151]). The i](twa) parameters are given in 15]. If we go on and calculate the 



diagonal matrix elements using the same procedure, we get an inadmissible large value for 
the gap. If we try the tight-binding parameters proposed by Papaconstantopoulos |16fl for 
Cu metal, we do not get the right gap as well. Also the Cu on-site parameters that reproduce 
correctly the electronic band structure of the superconducting perovskite. YBa 2 Cu 3 7 fail. 
Cu orbitals have an important influence on the gap edges in the electronic band structure 
of the Cu-based chalcopyrites. 

In the semiconducting Cu-based chalcopyrites, the s-like orbital plays a mayor role in 
fixing the lower edge of the conduction band while the p-like one influences the position 
of upper edge of the valence band. The d-like Cu-orbital mostly fixes the value of the 



chalcopyrite gap. These orbitals and the p-like Se ones repe 



valence band edge upwards so that the gap is diminished 
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each other and push the upper 



18, 



19j |. Consequently, we have 



fixed the Cu on-site parameters for the whole series in such a way that we get the lowest 
possible deviation from the semiconducting gap. More exactly, we have selected the three 

— E^^) 2 as a function of EL 



Cu on-site parameters so that ^ series (-E' series — Ug Series )" as a lunction ol £/ ser ies 1S minimum. 
, /MTI<5S are the experimental values of the gap. Small further adjustments of the anion (Se) 



E E 



p on-site parameter for each chalcopyrite allowed us to get the right experimental gap for 
the whole series. The experimental values that we used are quoted in the next Table I. 



Table I- Experimental optical gap 



Chalcopyrite 


l eV \ 


CuAlS 2 


3.49 


CuAlSe 2 


2.67 


CuAlTe 2 


2.06 


CuGaS 2 


2.43 


CuGaSe 2 


1.68 


CuGaTe 2 


1.23 


CuInS 2 


1.53 


CuInSe 2 


1.04 


CuInTe 2 


1.02 



or the whole series of Cu-based chalcopyrites considered 



to set the tight-binding parameters [1 



The on-site tight-binding parameters that we 
the ones obtained from the Harrison's formulas 



get in this way are compared in Table II to 



15] 



Parameter 


Harrison 


This work 


E s [eV] 


-6.92 


-14.55 


E p [e^] 


-1.83 


-2.22 


E d [eV] 


-20.14 


-16.97 


r d [i] 


0.67 


1.15 



Table II- The Cu on-site tight-binding parameters. The parameter is defined in [IJ 



III. THE HAMILTON! AN 



The Hamiltonian is labelled with the atom numbers in the following way. We take into 
account first nearest neighbors interactions only. The Hamiltonian matrix takes the form 

4 
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-"67 
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Se2 


#12 


#23 


#25 


#27 


#22 











Se4 


#14 


#34 


#45 


#47 





#44 








Se6 


#16 


#36 


#56 


#67 








#66 





Se8 


#18 


#38 


#58 


#78 











#88 



(4) 



The diagonal sub-matrices are 9x9 for Cu and 4x4 for In and Se. The Hamiltonian matrix 
is altogether 42x42. 

Obviously, H 33 = H 55 . These refer to Cu. Hn = H 77 , they refer to In. H 2 2 = H44 
= H 66 = H 88 which describe the Se atoms. The non-diagonal sub-matrices describe the 
first-nearest neighbors interactions. Their tight-binding parameters were computed from 
the Harrison's formulas The anion p-on-site parameter was adjusted further to get 

the exact experimental value (see below) With these data, the Hamiltonian can be built up 
straightforwardly 

IV. RESULTS 

A. The ideal case 

To get the experimental value for the gap (see Table II) in our calculated band structure 
we did a small (about 8%) further adjustment to the p-on-site parameter for the Se atom. 
The Harrison formula gives 9.53 eV to be compared with our 8.789 eV. The electronic band 
structure is presented in the next Fig. 1. There are 2 Cu, 2 In and 4 Se atoms in the unit 
cell. Each Se atom contributes with 2, each In with 3 and each Cu with ls+5d occupied 
electronic states and therefore 26 bands in the valence band. The rest of the 42 bands (16) 
appear as empty conduction bands. The band structure appears in the next Fig. 1. 

FIGURE 1 
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Both the top valence band and the bottom conduction band are approximately parabolic 
at T and therefore in some calculations the effective band approximation should be a good 
one. The semiconducting optical gap is direct and is calculated as the difference between the 
energies T lc and TfJ which is equal to 1.04 eV (a fitted value to the experimental one). The 
matrix element for dipolar transitions between these states, <T4 i ;|r|r lc > is different from 
zero along the z-axis (since z is an element of the T 4 representation; it is proportional to 
<r2 } |r 4 |r lc > and r 4 ®r! = r 4 ) but equal to zero at T (since <r£ ) |z|r£ ) >^<r 4 |r 4 |r 4 > 

and r 4 <g> r 4 = r\). 

1. The valence band 

Immediately below the singlet state T\J (at the top of the valence band), we find a doublet 
T^. The dipolar moment along z is different from zero (<r 5 ^|2;|r5 iy > ==>- < r 5 |r 4 |r 5 > and 
r 4 ® T 5 = r 5 ). Therefore the dipolar moment at the top of the valence band (a triplet in 
the zincblende parent compound) breaks into a zero dipolar moment at the top vfj and 

(2) 

a non-zero one at the doublet T^J . The operator representing the quadrupole moment is 
proportional to 3z 2 — r 2 which belongs to r\ and the products of the type <r x |ri|r x > are 
always different from zero since I\ ® T x = T x and so a non-zero quadrupole moment will 
exist for all the valence band states. 

The 26 bands that conform the valence band are grouped together into three sub-bands 
separated by two in-band gaps. The first one (A in Fig.l) separates the upper valence band 
(UVB) from the middle valence band (MVB) and the in-band gap B separates this band 
from the lower valence band (LVB), At the top of the UVB there is a singlet, r 4l/ , separated 

(2) 

from a doublet, T 5 J , by a crystal field splitting, A c j s , of 16 meV which is zero in the parent 
zincblende compound as we mentioned above. Notice that the doublet remains such from 
r — Z but splits from T — X. 

The chalcopyrite crystal field brakes the zincblende symmetry in several ways. First, 
there are two different cations instead of one which transforms the symmetry from cubic 
to tetragonal. Secondly, the anion can be found displaced along the center line of the 
tetrahedron that it forms together with the two different cations. But also the tetragonal 
symmetry (j- — 1 where c and a are the lattice parameters) is broken. We will deal with 
the effect of these distortions below. 
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a. The upper valence band (UVB) The splitting of the triplet in the zincblendes into 
the singlet rfj and the doublet vf} here, can be traced easily to be due to the presence a the 
two cations. This can be done with the tight-binding program replacing the parameters for 
Cu with those of In to get the bands for InSe or reversing the way we replace the parameters, 
we can get CuSe. In both cases the bands show a triplet on the top of the valence band at 
T.When the bands of the zincblende are compared to the ones of the chalcopyrite, we realize 
that two further splitting occur at the top of the valence band, one at Z (AZ in Fig. 1) and 
another one at X (AX). 

There are 10 bands in the UVB that lie from to -5 eV (the origin is set at the top of the 
valence band in T as it is customary). The main contribution comes from p-like Se orbitals. 
The details of the composition are in Fig.2 where the density of states (DOS) is shown. The 
shadow areas are proportional to the contribution of the orbital identified at the upper right 
corner. 

b. The middle valence band (MVB) The inner-band gap A is small (16 meV, see 
Figs. 1 and 2). The MVB contains 12 bands; 10 of them are d-Cu orbital contributions. 
The deepest band of this group runs from Z± v + Z^ v — * T^J — > X$ as shown in Fig.l. 

c. The lowest valence band (LVB) The deepest group of bands, the LVB, is separated 
from the MVB by a large gap from T^J to T 3u of about 4eV. The main contribution comes 
from the s-Se orbitals (see Fig.2 for more details). The upper band of this group is a singlet 
T 3u followed very closely by a doublet T^J . In the zincblende parent compound these bands 
are degenerate. This splitting is due to the presence of a second cation. The upper band of 
this group Z lv + Z 2v — > T 3v is doubly degenerate from Z — T but splits from T — X. 

FIGURE 2 

2. The conduction band 

The conduction band (CB) minimum runs from Z\ c + Z 2c — ► Ti c — ► which is a 
singlet all along ZTX. At Z, nevertheless the band is degenerate and splits to a higher-in- 
energy band that runs from Z\ c + Z 2c — ► r 3c — > X^ . At X the band is again degenerate 
but at r 3c it is a singlet. From Fig.2, we see that the CB is divided into two clearly defined 
sub-bands separated by an in-band gap of about 0.2 eV. Each sub-band presents two peaks. 
The DOS is quite higher in the upper part the spectrum. The lowest conduction band (LCB) 
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goes from roughly 1-5 eV. The upper conduction band (UCB) runs from about 5.2-12 eV. 
The lower peak of the LCB is composed mainly from s-Cu and s-Se orbitals in the 1-3.7 eV 
energy region and from s-In and s-Se in the higher energy region. In the low energy region 
of the UCB, the main contribution is from p-Cu orbitals and in the higher energy peak it is 
from p-In ones (see Fig. 2 for more details). 

3. Comparison with other work 

Jaffe and Zunger (JZ) made an ab initio calculation of the electronic band structure 
of the chalcopyrite CuInSe2- Sometimes ab initio calculations do not get the semiconducting 
gap right. JZ got 0.98 eV. The experimental value is 1.04 eV. But otherwise, our calculation 
compares well with theirs. The valence band width (T^J — T„) is 14.8 eV in our work and 
13.8 eV in theirs. Both calculations agree in the fact that the top of the VB which is a 
triplet in the zincblende parent crystal structure is split apart into a singlet at the top and 
a deeper doublet. The A c f s = —0.08 eV in JZ and -0.016 in this work. We will show below 
that the difference shrinks as an effect of the distortions (both anion and tetragonal which 

has calculated this crystal field splitting A c f s = 
for the ideal case and A c fs = +0.01 eV when the tetragonal distortion (- = 2.008) is taken 
into account. In his case the bands are in a reverted order which means that the doublet is 
on the top of the valence band. This can be associated to his neglect of the anion distortion. 
There are some differences in both calculations. For example, for the width of the UVB we 
get 5 eV while JZ get 4 eV. This is actually the origin for the difference in the overall VB 
width. It is worth mentioning that the in-band gap A differs substantially in both works 
even qualitatively. The two bands that define this gap are reverted in JZ's work giving a 
-0.01 eV while we get a broad in-band gap of 1.67 eV. The overall width of the MVB does 
not differ very much. We get 1.5 eV and JZ 1.37 eV. On the other hand, the in-band gap B 
is higher in JZ (7.39 eV) while we get 4.79. Some of the differences are shown in detail in 
Table V below. 
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B. Influence of distortions 



It os obvious that a correct description of the chalcopyrites must include the effect of 
distortions (both, the anion and the tetragonal one). In this section we will show the effect 
of this inclusion. The tetragonal distortion means that the ratio j- = 77 7^ 1 and the anion 
distortion means that the anion is not located exactly at a(|, |, j) and equivalent positions 
but rather at a(u, |, ?). The new positions are shown in the next Table IV. 



Nb. 


atom 


ideal 


distorted 


I 


In 


(0,0,0) 


(0,0,0) 


2 


Se 


(I 1 1) 


Ki?) fl 


3 


Cu 


(io,|) 


(|,0, f)a 


4 


Se 


(3 1 3N 
V4' 4' 8' 


(4' 2 M 5 4 ) a 


5 


Cu 


(0,0, 1) 


(|,|,0)a 


6 


Se 


(I I 

\4> 4' 8/ 


(l-«,|,2)a 


7 


In 


( l 2) 


(0,|,=)a 


8 


Se 


U' 4' 8^ 


(i 



Table IV Positions of the atoms in the unit cell in the ideal and distorted case 



We use the same Hamiltonian described above with the same basis on each atom and 
introduce the changes that describe the corresponding distortions. We can use the same 
symmetry group to construct the Hamiltonian since X-ray experiments indicate that the 



space group is preserved 



3 y, 0, y , 



261 ] in spite of the distortions present. In this 



case the Se-atom (the anion) will have its first nearest neighbors (the cations) at different 
distances. Namely, 



R 



In >5e 



R 



.12 



a\\u A + 



16 



Rse 



*Cu 



R 



2.3 



1 (1 + 77)2 

a\ (u Y + — 

1 v 2 J 16 



(5) 
(6) 



and the direction cosines that intervene in the calculation are given by: 



Vector 



1 



R 
R 



12 



2:-! 



ua 
gig 

R23 



m 



4Hi2 



4i?23 
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qa 
va 
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Tabla IV The direction cosines once distortion is taken into account. These are equal 

always to 4= in the ideal case. 

The phase factors are to be corrected accordingly. Now it is straightforward to construct 
the tight-binding Hamiltonian that includes both, the tetragonal as well as the anion dis- 
tortion. We present in the next Fig. 3 our results for the influence that the inclusion of 
distortion has on the electronic band structure of CuInSe 2 . 

FIGURE 3 

The general picture of the bands is quite similar. The VB present the same three groups 
separated by two in-band gaps. We do not obtain an important reduction of the total band 
width (13.79 eV to be compare to the ideal case value of 13.8 eV). The conduction band is 
formed again by two groups separated by a small in-band gap. Also as a general trend, we 
can say that in several aspects the result including the distortions agrees better with the JZ 
calculation. 

The first think to notice is that since the space group remains the same, the degeneracies 
at the high symmetry points Z, T, and X also remain as it is to be expected. On the 
other hand, we see that the distortions have a mayor or minor effect on almost all the 
bands. Distortion diminishes the optical gap and as a consequence we get 0.95 eV in closer 
agreement with the 0.98 eV of JZ . The experimental value to which we have fitted the 
on-site parameters in the ideal case is 1.04 eV. We could retouch the p-on site Se parameter 
to get the experimental gap value back. We present our results here without this further 
adjustment to leave sharp the picture of the effect of distortion. 

In the UVB, the crystal field splitting A cfs = -0.09eV (JZ got -0.08eV) in contrast 
to the ideal case value — 0.016eV. So the distortions have a very important influence in 
this crystal field splitting. If we take into account the result by Yoodee et al. |23( quoted 
above, then we can conclude that the most important factor that influences the crystal field 
splitting is the anion distortion. The AZ and AX splitting enhance due to distortions as 
can be seen from the enhancement of the top valence band details in Fig. 3. This sub-band 
broadens and as a consequence the in-band gap A shrinks. 

The next MVB appears to be still a very narrow band from the d-Cu orbitals mainly. The 
effect of distortions is here a minor one and it is related essentially to the slight shrinking 
of the total sub-band width which enhances but slightly the in-band gap B. Notice that the 
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band Z± v + Z*> v — ► F\J — > X \J becomes almost flat as a result of the distortions. 

In the LVB which contains 4 bands associated with s-Se orbitals the main effect is to 
invert the sign of the crystal field splitting which reverts the order between T 3u and T^J so 
that now the doubly degenerate band is on top. We summarize the effect of distortion in 
Table VI where we compare the energy values for several bands that we get for the ideal 
case and for the result after distortions are included with the ones obtained by JZ. 



State 


JZ 


Ideal 


Distortioned 


UVB-maximum 


(eV) 


(eV) 


(eV) 


F (2) 
1 Av 


0..00 


0.00 


0.00 


F (2) 


-0.08 


-0.02 


-0,09 


Z30 + z^ v 


-0.79 


-1.93 


-1.82 




-0.54 


-1.83 


-1.22 


UVB- minimum 


r (i) 

1 Av 


-5.15 


-8.11 


-7.89 


z± v + z$ v 


-5.12 


-7.98 


-7.87 




-5.13 


-7.92 


-7.82 


s-Se band 








r (i) 


-13.03 


-12.91 


-12.82 


^3v 


-13.06 


-12.90 


-12.90 


F ( 1 ) 
1 lv 


-13.83 


-14.80 


-14.79 


Z\ v + Z2V 


-13.00 


-12.91 


-12.82 




-13.46 


-14.00 


-13.99 




-13.20 


-13.36 


-13.38 




-13.31 


-13.57 


-13.47 


Other values 


Width band s-Se at Y 


0.80 


1.90 


1.96 


Gap A 


-0.01 


1.67 


1.23 


Gap B 


7.39 


4.79 


4.94 


Az 


0.5 


0.11 


0.27 


Ax 


0.4 


0.26 


0.37 




0.03 


-0.01 


0.08 
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Table V Comparation between energy values at some high simmetry points in the ideal 

and distorted cases. 

In the conduction band the major effect of distortion is to shrink the sub-band width of 
the s-Cu orbitals and to enhance the one corresponding to s-In orbitals. The in-band gap 
between the two enhances. We compare some CB values in the ideal case with the ones when 
the distortions are included in the calculation to show the effect of distortion and quote the 
corresponding values obtained by JZ in the next Table VI. 



State 


JZ 


Ideal 


Distorted 


Conduction band 


(eV) 


(eV) 


(eV) 


Tic = E g (gap) 


0.98 


1.04 


0.95 


r 3c 


3.24 


3.77 


3.25 




2.76 


3.14 


2.83 




2.25 


1.87 


1.56 



Table VI The effect of distortion in the Conduction Band. 
C. The (112) surface 

In the next Fig.4, we present our result for the effect that the distortions have on the 
electronic band structure of the (112) surface of CuInSe 2 . The open circles are the ideal case 
and the black ones are the result when distortions are included. In general, distortions do 
affect the VB and the CB in a noticeable way. 

FIGURE 4 

Of great interests the effect of distortions on the surface states. The full line in Fig.4 
represents the border of the bulk valence band (VB) and of the conduction band (CB). 
Three surface states denoted as 1SE, 2SE and 3SE in Fig.4 appear. The 1SE state starts 
at T somewhere in the middle of the semiconducting gap. Both 2SE and 3SE appear to be 
closer to the CB minimum. Somewhere in the interval K-J the 2SE state crosses the state 
1SE and becomes the one nearest to the top of the VB. If we follow 1SE further on, we 
see an important influence of distortions on this state. Indeed, in the ideal case, 1SE enters 
the CB region and becomes a resonance while if distortions are included, the state remains 
within the energy gap region as a pure surface state. A similar effect occurs with the 3SE 
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FIG. 1: Influence of distortions in the electronic band structure of the (112) surface of CuInSe2- 
Both, the tetragonal and the anion distortion are included. The open circles are the ideal result 
and the black ones the one with distortions included. 

state around K which is a resonance in the ideal case but becomes a pure surface state 
as an influence of distortions. The 2SE state is switched towards the top of the VB when 
distortions are included. So the (112) CuTnSe 2 surface shows three pure surface states as a 
consequence of distortions that would partially become resonances would the lattice remain 
ideal. We conclude that distortions have an effect also in the surface states that is important 
and cannot be neglected. 

V. CONCLUSIONS 

We have studied the effect of distortion in the electronic band structure of the chalcopyrite 
CuInSe2- We find that its effect is important and that the inclusion of both the tetragonal 
and the ionic distortion are important to get a proper description of the electronic bands. 



13 



Once these are included we find that that the tight-binding Hamiltonian gives an accurate 
enough result to be useful for further calculations of surfaces, monolayers and interfaces 
and more complicated system that include this material. We have used the Hamiltonian 
together with the Green's Function Matching Method to calculate the effect of distortion on 
the (112) surface of this material. We find that distortions have an important effect on the 
three surface states that we found on this surface which is used in important technological 
applications. 
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FIGURE CAPTIONS 

Fig.l The electronic band structure of CuInSe 2 . 

Fig.2 Contribution to the DOS from the different orbitals at different energies. 

Fig. 3 Influence of distortion on the electronic band structure of CuInSe2- We used u = 
0.224 and rj = 1.004 (the ideal values are 0.25 and 1 respectively). The full line is the ideal 
case and the dot lines represent the effect of distortion. 

Fig. 4 Influence of distortions in the electronic band structure of the (112) surface of 
CuInSe 2 . Both, the tetragonal and the anion distortion are included. The open circles are 
the ideal result and the black ones the one with distortions included. 
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